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ABSTRACT 


In  contradistinction  to  the  known  theory  on  complex  splines 
which  are  defined  on  the  boundary  of  a  region  ^in  <jT),  we  define 
complex  planar  splines  on  a  region  itself  as  a  complex  valued 
continuous  function  which  is  defined  piece-wise  on  suitable  meshes 
of  that  region.  The  main  idea  is  to  use  non-holomorphic  functions 
as  pieces,  since  holomorphic  pieces  would  lead  to  just  one  holo- 
morphic  function  on  the  whole  region  by  a  well  known  identity 
theorem  in  the  theory  of  functions  in  one  complex  variable. 
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Some  of  the  techniques  used  are  available  from  the  theory 
of  finite  elements.  But  also  are  considered  new  aspects,  namely 
mapping  properties  of  a  complex  planar  spline  v  and  the  differ¬ 
ence  f  -  v  where  f  is  in  general  a  holomorphic  function.  For 
triangular,  rectangular,  parallelogram  meshes  and  meshes  on  cir¬ 
cular  sectors,  explicit  expressions  are  provided  and  also  proper¬ 
ties  of  the  newly  introduced  complex  planar  splines  are  studied. 


1.  INTRODUCTION 

In  this  paper  we  are  concerned  with  the  approximation  of 
complex  valued  functions  by  functions  which  we  would  like  to  call 
complex  plana*  splines. 

In  the  current  literature  (see  e.g.  Ahlberg  fl],  Ahlberg- 
Nilson-Walsh  [2,3,4],  Atteia  [7],  and  Schoenberg  [12,13]),  complex 
splines  are  defined  on  the  boundary  of  a  given  region  and  are  then 
extended  into  the  interior  by  Cauchy's  integral  formula.  However, 
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this  extension  process  is  not  easy  to  execute  numerically.  There¬ 
fore,  we  offer  another  approach  here  which  in  spirit  originates 
from  the  theory  of  finite  elements.  (For  a  comprehensive  infor¬ 
mation  on  this  subject,  see  the  references  in  Schwarz  [15]) .  We 
subdivide  a  given  region  into  me*  he*  and  define  a  complex  valued 
function  on  that  region  piecewise  on  each  mesh.  The  functions 
defined  on  each  individual  mesh  will  be  called  element*.  The 
vertices  of  the  occurring  meshes  will  be  called  gn.id  point*  and 
the  set  of  all  grid  points  a  g nid. 
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The  minimum  requirement  we  impose  is  the  continuity  of  that 
piecewise  defined  function.  Any  such  piecewise  defined  complex 
function  which  is  continuous  will  be  called  a  complex  plana K 
ipliue. 

The  continuity  already  has  a  very  serious  implication.  If 
we  would  try  to  define  a  complex  planar  spline  by  holomorphic 
elements  like  polynomials,  then  by  the  well  known  identity  theorem 
(e.g.  Diederich-Remmert  [9,  p.  132,  Theorem  60])  all  the  elements 
represent  just  one  holomorphic  function.  The  consequence  is  that 
it  makes  no  sense  to  work  with  holomorphic  elements.  Therefore,  we 
have  to  use  nonholomorphic  elements. 

Very  simple  nonholomorphic  functions  are  polynomials  in  the 
complex  variable  z  and  its  complex  conjugate  z  .  These  func¬ 
tions  have  the  form 

P(Z'*>  =  J-0  zJik  '  aJ*  '  *  •  (1-1) 


The  number 


3p  = 


max 


3 ,k=0 ,1, . . . ,n 


{j  +  k : a  .  .  f  0} 

j 


(1.2) 


will  be  called  the  degree  of  p  .  If  3p  =  1  we  shall  say  that 
p  is  liman,  if  3p  =  2  we  shall  say  that  p  is  quadratic , 
for  3p  =  3,4,5  we  shall  use  the  words  cubic,  quantic,  quintic 
respectively . 
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If  we  use  the  representation 


z=x+iy,  z  =  x  -  iy 


(1.3) 


then 


x  =  i(z  +  z)  ;  y  =  jt-  (z  -  z) 


(1.4) 


That  means  a  function  in  the  real  variables  x  and  y  can  be 
transformed  into  a  function  depending  on  the  complex  variables 
z  and  z  ,  and  vice  versa. 

Thus,  a  function  u  =  u(z,z)  may  also  be  regarded  as  a 
function  in  x  and  y  .  If  u  is  continuously  dif f erentiabl 


:o  x 

and 

y  # 

then 

we  have 

u  ,  = 

u  + 

U-  ; 

u  = 

i  (u  -  u-)  , 

(1.5) 
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If  u  is  twice  continuously  differentiable,  then 
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(1.7) 


In  the  sequel  we  treat  triangular,  rectangular  and  paral¬ 
lelogram  meshes  and  meshes  on  circular  sectors.  The  elements 
defined  on  those  meshes  will  be  as  simple  as  possible.  Besides 
the  aspects  known  from  the  theory  of  finite  elements,  as,  for 
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instance,  interpolating  properties,  computational  aspects, 
and  an  error  analysis,  there  are  new  aspects  which  can  be 
summarized  under  the  word  mapping  pnopnntiai*  of  the  newly 
defined  complex  planar  splines.  These  new  aspects  concern 
the  following  questions,  among  others: 

(I)  How  close  are  complex  planar  splines  to  comformality? 

(II)  Are  complex  planar  splines  quasiconformal? 

(Ill)  Are  complex  planar  splines  open  mappings? 

(IV)  Is  the  boundary  maximum  principle  valid  for  complex 
planar  splines? 

In  some  cases  these  questions  apply  also  to  the  difference 
e  =  f  -  v  between  a  certain  function  f  and  a  complex  planar 
spline  v  . 

Since  the  interpolating  formulae  for  the  complex  case  look 
different  than  the  corresponding  formulae  for  the  real  case,  we 
believe  that  it  is  reasonable  to  state  these  formulae  here.  We 
shall  see  that  the  complex  interpolation  and  L^~ approximation 
problem  reduce  to  two  real  problems  such  that  no  new  error  analy¬ 
sis  is  needed.  But  the  situation  changes  for  instance  for  the 
uniform  approximation  problem. 

Besides  the  approximation  of  functions  by  complex  planar 
splines,  still  another  application  seems  possible,  namely  solving 
complex  differential  equations  like  Beltrami's  equation 

fs  *  ^  +  ^  ‘i-8’ 
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without  splitting  it  into  real  and  imaginary  parts.  The 
equation  (1.8)  was  recently  treated  numerically  by  Weisel 
r 161  for  the  case  y  =  0  by  solving  the  corresponding  real 
system  with  finite  element  techniques.  A  systematic  treat¬ 
ment  of  (1.8)  can  be  found  in  Wendland  [17].  Furthermore, 
for  the  conformal  mapping  problem  there  also  exists  an  invest¬ 
igation  by  Bosshard  [8]  on  the  use  of  finite  elements  again  by 
treating  corresponding  real  problems.  Also  in  this  case  a 
direct  attack  seems  possible  which  then  avoids  the  computation 
of  the  conjugate  function  to  the  computed  real  part  of  the  mapping. 

Still  another  application  is  the  automatic  construction  of 
meshes  with  desired  behavior.  In  many  cases,  for  instance,  it 
is  desirable  to  construct  meshes  which  concentrate  at  a  certain 
point.  An  example  is  shown  in  Figure  1,  where  the  range  of  an 
interpolating  complex  planar  spline  is  sketched. 


_ Figure  1 _ 

These  applications  provide  motivation  for  our  study  of  com¬ 
plex  planar  splines  and  their  basic  properties. 

2.  LINEAR  COMPLEX  PLANAR  SPLINES  ON  TRIANGULAR  MESHES 
A  triangle.  A  is  defined  by  its  three  vertices  P^,P2»P3 
which  are  supposed  to  be  three  pairwise  different  complex  numbers 


not  located  on  a  straight  line.  The  three  edges  of  A  will 


be  designated  by  ,  P2  P3  ,  P3  P^  where  the  order  of 

the  indices  is  irrelevant. 

We  first  investigate  one  single  element  on  a  triangle  A  . 

We  will  use  a  linear  element  of  the  form 

p(z,z)  =  a  +  bz  +  cz  ;  a,b,c  e  (f  .  (2.1) 

For  simplicity  of  the  notation  we  shall  write  p(z)  instead 
of  p(z,z)  . 

It  is  clear  that  the  interpolation  problem  p(P^)  =  Cj  , 
j  =  1,2,3  has  a  unique  solution  for  any  three  complex  numbers 
^l'^2'^3  *  This  solution  is  given  by 


a  >  <C1<P2P3 

-  W 

+  ?2(P1P3 

-  PjPj)  +  ^3(P1P^  -  P^P2))/6  , 

(2.2) 

b  =  U-l  (P2 

+  c 

2(P3  " 

+  C3(Pj  -  Pj))/6  , 

(2.3) 

c  =  (C1(P3  " 

P2)  +  i 

2(P1  "  P3) 

+  C3(P2  -  P1))/6  , 

(2.4) 

where 

6  -  21 

Im  (PjP^  + 

P2*I  +  P3PT)  • 

(2.5) 

(Note  that  Re(z)  and  Im(z)  mean  the  real  and  imaginary  parts 
of  z  respectively.) 

It  is  easy  to  see  that  any  three  points  z^,z2,z3  6  $  form 
a  triangle  if  and  only  if 
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Im{z1z2  +  z2z3  +  z-^z.^)  7*  0  (2.6) 

and  this  property  is  invariant  under  translation. 

For  the  btandand  tfidangle.  P^  =  0  ,  P2  =  1  ,  =  i  the 

above  formulae  reduce  to 


with 


a  =  r,1  , 

b  =  |  Ux<i  -  1)  +  C2  -  iC 3)  , 

c  -  j  (~C L  (i  +  1)  +  C2  +  iC3)  » 

6  =  -2i  . 


(2.7) 

(2.8) 

(2.9) 

(2.10) 


For  any  given  triangle  with  vertices  Pi'P2'P3  we  can 
construct  three  special  elements  Pj  ,  j  =  1,2,3  by  solving 
Pj  (pk)  =  ,  j,k  =  1,2,3  ,  where  is  the  common  Kronecker 

symbol.  These  elements  are  usually  called  jjo/tm  efemen-tA.  They 
have  the  property  that  the  general  interpolating  element  p  with 
p(Pj)  =  Cj  ,  j  =  1,2,3  can  be  written  in  the  simple  form 

3 

P  =  Sj  Pj  *  (2.11) 

Their  importance  lies  in  the  fact  that  they  can  be  used  to 
construct  a  basis  for  the  linear  space  of  all  complex  planar 
splines,  as  will  be  seen  later. 


Definition  2.1.  A  linear  element  p  =  a  +  bz  +  cz  for  any 
a,b,c  e  is  called  degenerate  if  |b|  =  |c|  . 

Let  R  be  a  region  in  £  and  f  :  R  -*■  a  mapping 
which  has  continuous  partial  derivatives  f  and  f—  .  If 
f  is  quasiconformal ,  the  numbers 


d(z)  = 


Vi  (z)  = 


fz(z) I  +  I fz(z) 

f  (z) |  -  |f-(z) 
z  '  1  z 


fi(z) 

f2(z) 


(2.12) 


(2.13) 


are  called  dilatation  quotient  and  complex  dilatation  at  the 
point  z  respectively  (Lehto-Virtanen  [10,  p.  52  &  191]).  We 
always  have  d(z)  >  1  .  In  case  d(z)  =  1  ,  z  could  be  called 
a  conformal  point.  If  d(z)  =1  for  all  z  ,  then  f  is  con¬ 


formal  . 


The  following  theorem  gives  some  properties  of  a  linear 


element  p 


Theorem  2.1.  A  linear  element  p(z)  =  a  +  bz  +  cz  ,  a,b,c  e  <T 
has  the  following  properties: 


1.  If  zi'z2  e  $  * 


then 


P  (  ( 1  “  +  \z2)  =  (1  -  X)p(z^)  +  Xp(z2) 


for  all  X  e  ]R 


(2.14) 


if  and  only  if  p  is  non-degenerate. 

3 .  p  is  an  orientation  preserving  homeomorphism  if  and 
only  if  J  b |  >  |  c|  . 

4 .  p  is  an  open  mapping  if  and  only  if  p  is  nondegenerate. 

5.  p  is  quasiconformal  if  and  only  if  |b|  >  |c|  .  In  that 
case  its  dilatation  is  a  constant  given  by 


d 


b 

b 


(2.15) 


and  the  complex  dilatation  is  a  constant  given  by 


y 


(2.16) 


6.  If  we  apply  p  to  an  angle  t  ,  0  <  t  <  n  f  the  angle 
t  is  distorted  by  the  angle 


t  =  arctan 


+  Re (be)  -  Im(bc)  tan  x 


Im  (be)  +  i( 


-  ] b |  ) tan  x 


1 , ,  ,2 

-  ^jb  +  c|  cot  x 


(2, 


Proof :  Properties  1-5  are  immediate.  The  proof  of  property  6 

follows  by  going  through  standard  computations. 

We  now  come  to  another  important  application  of  Theorem  2. 


Definition  2.2.  Two  triangles  A2 

m-Cqhbo filng  if  they  share  a  common  edge. 
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are  called  neighbor*  or 


liTsI»V«C:aul 


common  edge  P^'  P2  .  (see  Figure  2)  .  Further  let 


be  a  linear  element  in  A.  ,  j  =  1,2  .  If 
-  3  — 


p.  (P .)  =  p  (P  , )  ,  j  =  1,2 
J  ^  J 


(2.18) 


then 


p^(z)  =  p2(z)  for  all  z  e  P^  ,  and  consequently 


p(z)  = 


p^(z)  for  z  e  A^ 
p2(z)  for  z  e  A2 


(2.19) 


is  continuous  on  A^  u  A2  . 


Proof .  Follows  directly  from  Theorem  2.1  (property  1) 


Now  if  we  subdivide  a  region  into  triangular  meshes  and 
define  a  linear  element  in  each  mesh  we  obtain  a  complex  planar 
spline  if  we  impose  condition  (2.18)  for  each  pair  of  neighbor¬ 
ing  triangles.  In  the  triangulation,  however,  we  do  not  allow 
that  a  vertex  of  any  triangle  is  interior  to  any  other  edge. 
More  specifically  the  triangulation  has  to  be  ptopct  (Prenter 
111,  p.  127  1)  . 
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The  form  elements  introduced  earlier  are  used  to  construct 
a  basis  for  the  linear  space  V  of  all  complex  planar  splines. 
Let  R  be  a  region  in  ^  ,  subdivided  into  finitely  many  tri¬ 
angles  and  P  ,  P2,...,  its  grid  points.  A  complex  planar 

spline  associated  with  the  gridpoint  and  defined  by 

1  for  z  =  P.  , 

0  for  z  e  { P  ^ ,  P2,  ...,  P^l  -  (P.l  ,  j  =  1,...,N 

will  be  called  a  ( global )  faosim  function .  It  can  be  constructed 
piecewise  from  the  (local)  form  elements  already  known.  Assume 
that  T. , T„ , . . . , T,  are  the  triangles  which  have  P.  as  a  common 

1  Z  Kj  I 

vertex.  Then  combine  those  form  elements  defined  in 

Tf , T2 , . . . , which  have  the  value  one  at  P^  and  zero  at  the 

two  other  vertices  with  the  zero  elements  on  all  other  triangles. 

Clearly  V  =  <vl ,  v2  - - -  vN>  . 

A  complex  planar  spline  in  general  does  not  represent  an 
open  mapping  even  if  all  elements  are  not  degenerate.  To  see 
this  we  define  a  complex  planar  spline  on  the  square 

Q={z:0^Rez5l,  0  s  Im  z  <  1}  ,  (2.20) 

which  we  divide  into  two  triangles 

A  ^  =  { z  f  Q  :  Rez  +  Imz*'!},  (2.21) 


{ z  c  Q  : 


Re  z  +  Im  z  >  1 }  , 


(2.22) 
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ii 


and  define 


f  z 


for  z  e  A^  , 


v  =  < 


(2.23) 


J  1  +  i  -  iz  for  z  e  A2 

which  is  then  a  complex  planar  spline  on  Q  . 

However,  the  range  of  Q  under  v  is  such  that  the 

range  of  any  open  set  in  Q  containing  parts  of  the  diagonal 
of  Q  is  not  open. 

But  we  have  the  following  Lemma  (where  proof  is  obvious) . 


Lemma  2.1.  If  a  complex  planar  spline  represents  a  univalent 
function  it  is  an  open  mapping. 

A  consequence  of  this  Lemma  is  the  following. 


Theorem  2.3.  Let  v  be  a  complex  planar  spline  which  inter¬ 
polates  a  univalent,  holomorphic  mapping  f  on  the  grid  points 
of  a  sufficiently  fine  triangular  grid  which  is  inside  of  the 
domain  of  definition  of  f  .  Then  v  is  univalent  and  open. 

Proof .  Since  f  is  conformal,  a  small  triangle  is  mapped  such 
that  the  images  of  the  three  vertices  form  a  triangle.  There¬ 
fore,  the  interpolating  complex  planar  spline  maps  the  triangular 
grid  onto  another  triangular  grid  in  a  univalent  way.  The  result 
follows  from  Lemma  2.1. 


Lemma  2.2.  Let  f  be  a  holomorphic  function  in  a  region  R  and 
g  a  holomorphic  function  in  R  =  {z  :  z  e  R)  .  Define  a  function 
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h  on  R  h(z)=f(z)+g(z),zeR.  If  h  is  not  con¬ 

stant,  then  | h |  does  not  admit  a  maximum  in  R  . 

2 

Proof :  It  is  sufficient  to  show  that  |h|  has  no  maximum  in 

R  .  If  A  represents  the  Laplace  operator  as  defined  in  (1.7) 
we  obtain 


2 

But  this  implies  that  |h|  is  subharmonic.  Since  h  is  not 
constant,  |hj  does  not  admit  a  maximum  in  R  (Ahlfors  [6, 
p.  245]) . 

Particularly  a  non-constant  linear  element  p  itself  has 
the  property  that  it  admits  no  maximum  in  the  interior  of  the 
triangle  in  which  it  is  defined.  This  is  a  stronger  property 
than  property  4  of  Theorem  2.1. 

Clearly,  a  complex  planar  spline  is  not  subharmonic  in  gen¬ 
eral.  As  an  example  take  an  interpolating  spline  which  is  1 
at  one  interior  gridpoint  and  zero  at  the  other  grid  points. 


3 .  QUADRATIC  COMPLEX  PLANAR  SPLINES  ON 
RECTANGULAR  AND  PARALLELOGRAM  MESHES 

First  we  study  splines  on  rectangles  whose  sides  are 
parallel  to  the  axes.  Such  a  rectangle  Q  (called  quabla 
in  Physics)  is  defined  by  its  four  vertices  Pi'P2,P3'P4  which 
are  to  be  understood  as  pairwise  distinct  complex  numbers  with 
Im  P^  =  Im  P2  ,  Im  P3  =  Im  P4  ,  Re  P^  =  Re  P^  ,  Re  P2  =  Re  P3 
arranged  in  positive  orientation. 
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We  first  investigate  one  single  element  on  a  rectangle 
A  quadratic  element  of  the  form 

_  2  _ _2 

p(z)  =  a  +  bz  +  cz  +  d(z  -  z  )  ,  a,b,c,d  e  <t  (3.1) 


will  be  used.  That  the  use  of  this  element  is  reasonable  will 
be  seen  in  Theorem  3.1.  The  interpolation  problem 
p(Pj)  =  Cj  ,  j  =  1,2,3, 4  again  has  a  unique  solution  for  any 
four  complex  numbers  ^  i 2 3 4  ,  and  this  is  as  follows: 


=  (1/6) ((P,  - 


P3}  C1  "  (P4 


+  (Pt  - 


P^)  C 
4'  ^2 


^  3 


<P2  *  P2>  ?4>  ' 


(3.2) 


b  =  (2/6)  (~P3  l,1  +  P4  C2  ”  P1  ^3  +  P2  '  (3.3) 

c  =  (2/6) (P3  ^  "  p4  +  P1  C3  ‘  P2  C4}  '  (3*4) 

d  =  (1/6)  (^  -  C2  +  ^3  “  ^4  )  '  (3.5) 

where 

"  4(P2  -  P1)(P3  '  P2)  • 
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(3.6) 
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In  the  case  of  quadratic  complex  planar  splines  on  rectan¬ 
gular  meshes,  we  have  the  following  analogue  of  Theorem  2.1. 

Theorem  3.1.  Let  P,Q  any  two  distinct  points  of  £  and 
p  a  quadratic  element  as  defined  in  (3.1)  with  d  ^  0  . 

Then 

p((l  -  X ) P  +  XQ)  =  (1  -  A)p(P)  +  Xp(Q)  for  all  X  £  [0,1]  (3.7) 

if  and  only  if  the  straight  line  through  P  and  Q  is  parallel 
to  the  x-  or  y-  axis . 

Proof .  Set  z^=  (1  -  X)P  +  XQ  .  Then  because  of  Theorem  2.1 
and  d  f-  0  ,  (3.7)  is  equivalent  to 

z*  -  z2  =  (1  -  X)  (P2  -  P2)  +  X  (Q2  -  Q2)  for  all  X  e  [0,1].  (3.8) 

Let  P  =  x  +  iy  and  Q  =  u  +  iv  .  After  routine  computations 
we  deduce  from  (3.8)  that  x(v  -  y)  =  u(v  -  y)  which  yields  the 
assertion . 


Definition  3.1.  Two  rectangles  will  be  called  neighbor*  or 
ne.igkbon.ing  if  they  share  a  common  edge  (Figure  3)  . 


Theorem  3.2.  Let  ] _ |^,  be  two  neighboring  rectangles 

whose  sides  are  parallel  to  the  axes,  with  the  common  edge 


P1P2  (Figure  3).  Let  be  a  quadratic  element  of  the  form 

(3.1)  defined  on  [ ,  j  =  1,2  .  If 

t] 


then 


pl(Pj)  =  p2  '  J  =  2 


(3.0) 


P1(z)  =  p2(z)  for  all  z  e 


Consequently 


p(z)  =  < 


Px(z)  for  z  e  PI  1 
P2 (z)  for  z  e  1  |  „ 


is  continuous  on  I  1^  u  1  t 


Proof:  follows  directly  from  Theorem  3.1. 


(3.10) 


Now  if  we  subdivide  a  region  into  rectangular  meshes  where 
the  sides  of  the  rectangles  have  to  be  parallel  to  the  axes 
then  we  obtain  a  complex  planar  spline  if  we  impose  condition 
(3.9)  for  each  pair  of  neighboring  rectangles. 


Theorem  3.3.  Let  R  be  any  region  in  <£  ,  fa  holomorphic 
function  in  R  and  p  a  quadratic  element  as  defined  in  (3.1). 
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If  f  -  p  is  not  constant,  then  j f  -  p[  does  not  admit  a 


maximum  in  R 


Proof:  Since  f  -  p  is  a  sum  of  a  holomorphic  function  in 

z  and  a  holomorphic  function  in  z  ,  Lemma  2.2  applies  and 
yields  the  desired  result. 

This  theorem  applies  particularly  to  p  itself. 

More  information  about  the  quadratic  element  p  as  defined 
in  (3.1)  can  be  deduced  from  its  Jacobi  determinant 


J  "  IPZI  2  -  IPil2 


(3.11) 


(see  Lehto-Virtanen  [10,  p.  136])  which  reads  here: 


J(z)  =  |  b  +  2dz|  2  -  |  c  -  2dz|  2 

=  |  b|  2  -  |  c(  2  +  2(db  +  dc)  z  +  2(db  +  dc)z  . 


(3.12) 


The  set 


H  =  {z  :  J(z)  =  0} 


(3.13) 


is  either  a  straight  line  or  a  point  in  <{:  . 


Definition  3.2.  The  quadratic  element  p(z)  =  a  +  bz  +  cz 
2  -  2 

+  d(z  -  z  )  ,  a,b,c,d,z  c  (f  is  called  de^em-ratc  if 


=  db  +  c(c  =  0 


(3.14) 


Now,  H  of  (3.13)  is  a  straight  line  in  <p  if  and  only  if  p 
is  not  degenerate.  If  p  is  degenerate  then  (3.14)  implies 
d  =  0  or  j b  j  =  | c  |  If  p  is  not  degenerate  then  the 
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straight  line  H  forms  the  angle  a  with  the  real  axis  which 
is  given  by 


tan  a  = 


Re  (db  +  dc) 
Im  (db  +  dc) 


0  <  a  <  tt 


(3.15) 


Let  us  assume  now  that  p  is  not  degenerate.  If  the  domain  of 
definition  of  p  is  any  compact  set  S  in  <£  which  is  located 
in  the  half  plane 

H+  =  {z:  J(z)  >  0}  (3.16) 


then  p  is  locally  an  orientation  preserving  homeomorphism  in 
S  which  then  is  also  quasiconf ormal  in  S  since  its  dilatation 
quotient  (see  2.12)  is  bounded. 

In  order  to  find  out  whether  p  is  a  global  homeomorphism 
we  study  the  solutions  of 

P ( z 2 )  -  p(z1)  =  0  .  (3.17) 

If  we  use  the  abbreviations 


x  =  z2  -  ,  y  =  z2  +  (3.18) 

equation  (3.17)  reads 

p { z 2 )  -  p(z^)  =  x(b  +  dy)  +  x(c  -  dy)  =  0  .  (3.19) 

From  this  it  follows  that 

d ( p ( z j )  -  p ( z 1 ) )  +  d (p ( z 2 )  -  p(z1))  =  xo  +  xo  =  0  ,  (3.20) 


where  o  was  already  introduced  in  (3.14).  If  we  use  this  equa¬ 
tion  to  eliminate  x  from  (3.19)  we  obtain 
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ya  +  y5  +  J  b ) 2  -  | c ) 2  =  0  (3.21) 

in  case  x  0  .  Let  J(z^)  >0  ,  j  =  1,2  .  Then  from  3.12 
by  forming  J(z^)  +  it  follows  that 

2 ( | b J  2  -  | c | 2 )  +  2oy  +  2oy  >  0  ,  (3.22) 

which  contradicts  (3.21). 

To  summarize:  If  p  is  not  degenerate  then  it  is 
univalent  in  both  half  planes  H+  and  H  =  {z  :  J(z)  <  0}  . 

Since  the  domain  of  definition  of  p  is  a  rectangle  R 
whose  sides  are  parallel  to  the  axes,  one  can  find  out  whether 
R  c  H+  just  by  inserting  a  suitable  vertex  into  J  .  To 
explain  that  we  distinguish  four  cases  according  to  the 
special  location  of  H+  in  $  . 


Case 

1  : 

7T 

2 

< 

a 

<  TT 

and 

H+ 

is  a  right  half  plane  in  (f 

Case 

_2 .  : 

a 

= 

0 

and 

H+ 

is 

an  upper  half  plane  in  ({:  , 

or 

0 

< 

a 

<! 

and 

H+ 

is  a  left  half  plane  in  (f  . 

Case 

_3  : 

TT 

2 

£ 

a 

<  TT 

and 

H+ 

is  a  left  half  plane  in  <f  . 

Case 

_4  : 

a 

= 

0 

and 

H+ 

is 

a  lower  half  plane  in  (f  , 

or 

0 

< 

a 

<  ! 

and 

H  + 

is  a  right  half  plane  in 

The  words  f  i>  f,t ,  light,  lo we* ,  uppzi  half  plane  are  used  in  the 
ordinary  sense. 
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Theorem  3.4.  Let  R  be  a  rectangle  whose  sides  are  parallel 
to  the  axes.  Call  p.^  the  lower  left,  P2  the  lower  right, 

P3  the  upper  right  and  P4  the  upper  left  vertex  of  R  . 

Further  let  p  be  a  non-degenerate  quadratic  element  ( as 
defined  in  (3.1))  on  R  .  Then  the  mapping  p  is  an  orienta¬ 
tion  preserving  homeomorphism  and  quasiconf ormal  on  R  if  and 
only  if  J(Pj)  >  0  where  j  is  determined  by  the  case  number 
j  to  which  H+  belongs,  j  e  (1,2, 3,4}  . 

Proof:  If  we  are  in  case  j  ,  j  e  {1,2, 3, 4}  ,  then  J(P^)  >  0 

is  equivalent  to  R  c  H+  . 

If  for  mnemotic  reasons  one  would  like  to  give  that  theorem 
a  name,  then  "  fiou/i  co/lne-t  th&otizm”  seems  to  be  very  suitable, 
since  all  four  corners  of  the  rectangle  R  are  involved. 

Under  the  assumptions  stated  p  will  be  an  orientation 
inv znting  homeomorphism  on  R  if  and  only  if  j(Pk)  <  0  where 
k  =  (j  +  2)mod  4  and  j  is  determined  as  before. 

If  a  rectangle  is  subdivided  in  this  way  into  m  •  n  little 
rectangles,  then  there  are  m  *  n  elements  having  altogether 
p  =  4  m  •  n  parameters.  Further  there  are  s  =  3mn-m-n-l 
continuity  conditions  leaving  p  -  s  =  (m  +  1) (n  +  1)  parameters 
free  where  (m  +  1) (n  +  1)  is  also  the  number  of  gridpoints. 

Now  we  can  adjust  the  element  (3.1)  to  the  case  where  the 
rectangle  has  any  position  in  the  plane.  If  one  of  the  edges  of 
a  given  rectangle  forms  the  angle  a  with  the  x-axis  then  instead 
of  (3.1)  one  has  to  use 


i 
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p(z)  =  a  +  bz  +  cz  +  d(z^  -  e^az^)  .  (3.23) 

Since  a  parallelogram  can  be  mapped  by  a  linear  transformation 
Z  of  the  type  (2.1)  onto  a  rectangle  whose  sides  are  parallel  to 
the  axes  one  can  also  work  with  parallelograms. 

If  ctjy  «2  e  with  ^  Aot2  for  any  X  e  IR  ,  describe  the 

directions  of  the  parallelogram  grid,  which  means  that  the  two 
angles  Bj  ,  j  55  1,2  with  the  real  axis  are  given  by 
tan  B,  =  Re  m^/Im  ou  ,  j  =  1,2  then  the  fore-mentioned  linear 
transformation  Z  yields  an  element  of  the  form 

p(z)  =  a  +  bz  +  cz  +  dfa^z  +  a^z)  (a2z  +  ct^z)  ,  a,b,c,d  e  (|)  .  (3.24) 

4.  SPLINES  ON  DISKS  AND  CIRCULAR  SECTORS 

Let  S  be  a  circular  sector.  If  we  divide  the  radius  into 
k  subintervals  and  the  opening  angle  into  Z  -  1  subintervals, 
we  obtain  in  total  k ( £  -  1)  meshes  of  two  types  which  we  would 
call  KZQ.tanguLi.an.  and  tx.ia.ng atari  meshes  respectively.  The  tri¬ 
angular  meshes  contain  the  origin  0  of  S  whereas  the  rectan¬ 
gular  meshes  do  not  contain  the  origin  (see  Figure  4). 


Figure  4 


It  is  straightforward  to  use  polar  coordinates  with  respect 
to  the  origin  0  of  S  in  this  situation.  Therefore  for  rec¬ 


tangular  meshes  we  use  an  element  of  the  form 

p(z)  =  pCre1^)  =  a  +  br  +  C0  +  dr<j>  ,  a,b,c,d  e  £  .  (4.1) 

Let  Pi'P2'P3'P4  be  the  vertices  of  one  specific  rectangular 
mesh  in  positive  orientation  such  that  $ ^  =  arg  P^  =  arg  P2  , 

<t2  =  arg  P3  =  arg  P4  ,  ^  =  |  P-L|  =  |  ?4  [  ,  r2  =  |P2|  =  |  P3| 

(see  Figure  5.) 

Then  the  solution  of  the  interpolation  problem 


p(Pj)  =  Cj  ,  e  <$  ,  j  =  1,2,3, 4 

(4.2) 

is  given  by 

a  = 

<r2«>2Cl  '  r1^2^ 2  +  rl*l*3  '  r2^4)/6  ' 

(4.3) 

b  = 

(_4,2C i  +  ^2^2  "  *1^3  +  ()>1^4)/(5  ' 

(4.4) 

c  = 

and 

(-r2Ci  +  rLC2  -  r1z3  +  r^^/6  , 

(4.5) 

d  = 

( r,  1  ~  C2  +  ^3  “  ^4  )/6  ' 

(4.6) 

i 

i 


where 
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_ Figure  5 _ 

Now  assume  that  0,P,Q  are  the  vertices  of  a  triangular 
mesh  in  positive  orientation  where  0  is  the  origin  of  the 
sector  S  (See  Figure  5) .  In  order  to  find  out  what  type  of 
element  to  use  here,  we  study  the  interpolation  problem  for 
a  rectangular  mesh  where  -*•  0  ,  -*•  0  and  =  ^4  '  such 

that  | P^ |  =  1 P^ 1  and  arg  P^  =  arg  P2  ,  arg  P^  =  arg  P^  . 

After  some  computation  we  find  c  0  in  (4.1).  This 
means  that  we  have  to  use  an  element  of  the  form 

p(z)  =  a  +  Br  +  yr<j>  ,  a,B,Y  e  .  (4.8) 

If  we  assume  that  P  and  Q  have  the  polar  coordinates 
(r^,<{>^)  and  (r2,<t>2)  respectively,  then  the  interpolation 
problem 

p(0)  =  ra  ,  p(P)  =  ,  p  (Q)  =  S3  ,  t,ir  t,  C3  e  <t  (4.9) 

has  the  solution 


"  ''l 

=  (-(*2  -  +  4^2  -  ♦1t3)/6 


(4.10) 


B 


(4.11) 
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and 


Y  =  (-Q2  +  ^3)/6 


(4.12) 


where 


5  =  l/r1(4>2  -  4^) 


(4.13) 


Theorem  4.1.  Let  S  be  a  circular  sector  subdivided  into 


meshes  as  described  above  and  m^  ,  j  =  1,2  two  neighboring 


me 


shes  sharing  the  common  edge  E  =  •  Let  be  an 


is  a  rectangular 

J 

mesh  or  it  has 

j 

the  form  (4.8)  if 

J 

m  .  is  a 

triangular  mesh. 

j  =  1,2  . 

3 

lf_ 

Pl(pj)  =  p2(P.) 

,  j  =  1,2 

(4.14) 

then 

p1(z)  =  p2(z) 

for  all  z  e  E  . 

(4.15) 

Consequently , 

px(z) 

p  (  z)  =1 

p2(z) 

l 

for  z  e  m^ 

for  z  e  m2 

(4.16) 

is  continuous  on 

ni  i  u  ^  / 

Proof:  The  elements ( 4 . 1)  on  : ectangular  meshes  as  well  as  the 
elements  (4.8)  on  triangular  meshes  are  linear  on  the  edges  of 
their  respective  domain  of  definition. 


The  element  p  introduced  in  (4.1)  has  the  same  form  as 
the  element  p  defined  in  (3.1)  and  used  for  ordinary  rectangles. 

In  order  to  see  this,  one  has  only  to  identify  Re  z  with 
r  and  Im  z  with  $  .  The  consequence  is  that  no  particular 
analysis  besides  the  analysis  for  elements  on  ordinary  rectangles 
is  required. 

We  end  this  section  with  the  computation  of  the  number  of 
free  parameters  in  a  complex  planar  spline  on  a  circular  sector. 

The  k  ( £  -  1)  meshes  distribute  in  (k  —  1)(£  —  1)  rec¬ 
tangular  and  l-l  triangular  meshes.  In  order  to  make  a  piece- 
wise  defined  function  continuous  we  have  to  impose  £(3k  -  1)  -  4k 
conditions.  The  number  of  parameters  is  (£  -  1) (4k  -  1)  leaving 
S,k  +  1  parameters  free  which  is  also  the  number  of  grid  points. 

5 .  LEAST  SQUARE  APPROXIMATION  AND  INTERPOLATION 
WITH  COMPLEX  PLANAR  SPLINES 

If  we  want  to  approximate  a  complex  valued  function  by  com¬ 
plex  planar  splines  of  a  certain  type,  we  end  up  with  minimizing 
a  real  valued  functional  defined  in  (£N  ,  N  c 

To  treat  such  a  problem  it  is  not  necessary  to  rewrite  it 
in  real  form.  Let  us  assume  that  we  have  to  handle  the  problem 

,N 

<t>(a)  =  min  ,  4>  :  (f  IR  (5.1) 

where  1  may  also  explicitly  depend  on  a  .  If  $  has  con¬ 


tinuous  partial  derivaties  with  respect  to  all  components  of 
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a  ,  then 


;  =  0  ,  j  =  1,2,..  .,N  ,  a  =  (alfa2,  .  .  .  ,an)  (5.2) 

is  a  necessary  condition  for  a  being  a  minimum  of  <p  .  This 
follows  immediately  from  (1.5)  and  (1.6) 

If  g,h  :  -*■  are  complex-valued  functions  possessing 

continuous  partial  derivatives  with  respect  to  a  and  a  , 
then 


(g)a  =  gi  '  (g) a  =  ga  (5-3) 

(gh)a  =  gah  +  ghl  (5.4) 

•^(g|2  -  ^  <g?'  -  gj  +  ggj  -  4  |g|2  ,  (5.5) 

and  if  g  does  not  depend  on  a  explicitly,  then 


9  a 


=  gag  • 


(5.6) 


The  least  square  problem  can  be  treated  along  the  lines 
of  Schultz  [14,  Ch.  6]  as  follows. 

Let  R  be  a  compact  set  in  (f  subdividable  in  meshes 
of  the  discussed  form,  f  e  L2(R)  and  V  the  linear  space  of 
all  complex  planar  splines  where  V  =  <v^ ,  v^  ,  vN>  and  the 

Vj  are  the  global  form  functions  defined  earlier  ,  j  =  1,2,...,N 
It  should  be  noticed  that  the  form  functions  are  by  definition 
real . 

The  problem  here  is  to  minimize 
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1  2  .N 

<j>  (a)  =  J|f(z)  -  l  a.v  .  |  dxdy  ,  a  =  (a^a^  .  .  .  ,  aN)  e  <f  (5.7) 

R  j=l  33 


Using  (5.2)  to  (5.7),  we  obtain 


N  _ 

<t>  (a)  =  J(  l  5\  v.  v  -  v  f ) dxdy  =  0 
R  k=l  K  3  K  3 


which  reduces  to  the  linear  system 


Ca  =  r  (5. C) 


where 


C  =  (cjk5  =  fvj  Vk  dxdy  »  3'  k  =  ,  (5.9) 

R 

r  =  (r.)  =  If  v.  dxdy  ,  j  =  1,2,. . . ,N  .  (5.10) 

J  R  J 

Since  C  is  a  real  matrix,  the  system  (5.8)  can  be  parti¬ 
tioned  into  the  two  real  systems, 

C  Re  a  =  Re  r  ,  C  Ima  =  Im  r  ,  (5.11) 

where  Re  a  is  the  vector  of  the  real  parts  of  a  ;  analogous 
meanings  apply  to  Re  r  ,  Im  a  ,  and  Im  r  . 

Clearly  (5.8)  has  a  unique  solution  since  is  strictly 

convex . 

Let  v  be  the  best  least  square  approximation  of  f  .  The 
error  analysis  can  be  directly  taken  from  the  real  case  (e.g. 
Schultz  [14,  Ch.  6])  since 
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II  f  -  v|!  =  1 1  Re ( f  -  v)  +  i  Im  (f  -  v)  II 

<  ||  Re  f  -  Re  v  ||  +  ||  Im  f  -  Im  v  || 


(5.12) 


This  means  that  the  order  of  convergence  is  the  same  as 
in  the  real  case,  but  the  convergence  constants  have  to  be 
doubled,  provided,  of  course,  that  Re  f  and  Im  f  are  of  the 
same  smoothness. 

An  interpolation  problem 


v(zk)  =  f ( z k )  ,  k  =  1,2, ...,K  (5.13) 

where  v  is  a  complex  planar  spline,  f  a  given  function,  and 
zk  ,  k  =  1, 2 , . . . , K  given  grid  points  may  be  as  well  partitioned 
into  two  real  problems  by  splitting  equation  (5.13)  into  real 
and  imaginary  parts.  The  above  remarks  on  the  error  analysis 
therefore  apply  also  here. 


6.  A  NUMERICAL  EXAMPLE 

We  subdivide  the  standard  triangle  0,1, i  in  the  usual  way 

Jr 

by  dividing  its  two  smaller  sides  into  1/h  =2  ,  k  =  0,1,..., 5 

pieces  of  equal  length  and  divided  the  resulting  little  squares 
of  side  length  h  diagonally  by  parallels  through  the  hypothenuse 
of  the  standard  triangle.  We  obtained  a  complex  planar  spline 
v^  by  interpolating  the  exp<.  rn  itial  function  at  the  grid  points 
of  that  triangle. 

In  the  following  Table  1  we  list  the  computed  values 


30. 


e  =  ||  exp  -  v,  ||  ,  the  corresponding  numerical  convergence 

n.  n  00 

order  c.  of  e,  ,  the  maximal  dilataion  quotient  d,  over 
h  h  n 

all  meshes  (compare  Formula  (2.15)),  and  the  corresponding 
convergence  order  c^  of  d^  . 


Table  1:  Spline  interpolating  the  exponential  function 
on  standard  triangle. 


h 

eh 

Ch 

dh  1 

5h 

1 

0.3730 

1.54 

2.229 

1.45 

1/2 

0.1285 

1.79 

1.4  51 

1.19 

1/4 

0.0372 

1.90 

1.198 

1.09 

1/8 

0.00996 

1.95 

1.093 

1.05 

1/16 

0.C0257 

1.98 

1.045 

1.03 

1/32 

0.000653 

1 

1.022 

The  fact  that  e^  approaches  zero  with  order  two  is 

of  course  known  (Schultz  [14,  Ch.  2]).  According  to  our 

computation  the  maximal  dilatation  approaches  one  with  order 

one.  The  number  d,  -  1  could  be  called  deviation  from  con- 

n 

formality.  To  the  best  of  our  knowledge,  neither 

numerical  values  for  the  dilatation  d,  nor  theoretical 

n 

investigations  on  the  behavior  of  d^  as  h  -*  0  exist  in 


literature. 
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